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We study the low-temperature thermal conductivity (re) of Dy2Ti20r and Yb 2 Ti 2 C >7 single crystals 
in magnetic fields up to 14 T along the [111], [100] and [110] directions. The main experimental 
findings for Dy 2 Ti 2 C >7 are: (i) the low-T k(H) isotherms exhibit not only the step-like decreases at 
the low-field (< 2 T) magnetic transitions but also obvious field dependencies in high fields (> 7 T); 

(ii) at T < 0.5 K, the k(H) curves show anisotropic irreversibility in low fields, that is, the k(H) 
hysteresis locates at the first-order transition with H || [100] and [110], while it locates between two 
successive transitions with H || [111]; (iii) the re in the hysteresis loops for H || [100] and [110] show 
an extremely slow relaxation with the time constant of ~ 1000 min. The main experimental findings 
for Yb 2 Ti 2 07 are: (i) the zero-field re(T) show a kink-like decrease at the first-order transition (~ 

200 mK) with decreasing temperature; (ii) the low -T k(H ) isotherms show a decrease in low field and 
a large enhancement in high fields; (iii) the low-T 1 re(IT) curves show a sharp minimum at 0.5 T for 
H || [110] and [111]. The roles of monopole excitations, field-induced transitions, spin fluctuations 
and magnetoelastic coupling are discussed. 

PACS numbers: 66.70.-f, 75.47.-m, 75.50.-y 


I. INTRODUCTION 

Rare-earth titanates R 2 Ti 2 07 (R = rare earth) have 
attracted extensive research interests due to their exotic 
magnetism.— These materials have a pyrochlore crystal 
structure with the space group Fd3m., in which the mag¬ 
netic rare-earth ions form a network of corner-sharing 
tetrahedra and are prone to a high degree of geometric 
frustration. However, these materials are greatly sensi¬ 
tive to weak perturbations (e.g. single-ion anisotropy, 
dipolar interaction or quantum fluctuations) beyond the 
nearest-neighboring exchange, which results in uncon¬ 
ventional low-temperature magnetic and thermodynamic 
properties. One famous phenomenon is the classical spin- 
ice state in Ho 2 Ti 2 0 7 and Dy 2 Ti 2 07 (with effective ferro¬ 
magnetic exchange and Ising anisotropy). 2-6 At low tem¬ 
peratures, their moments have an Ising anisotropy due 
to the strong crystal field with the local easy axis along 
the [111] axis. The ground states have a macroscopically 
degenerate “2-in, 2-out” spin configuration in each tetra¬ 
hedron. An interesting finding in the spin-ice materials 
is that the magnetic excitations could be emergent mag¬ 
netic monopoles— — Once the flipping of a spin occurs, 
a local “3-in, 1-out” or “3-out, 1-in” spin configuration 
forms, which is equivalent to yielding two opposite mag¬ 
netic monopoles in the adjacent tetrahedra— 


For Dy 2 Ti 2 07 , the external magnetic field along dif¬ 
ferent directions can easily break the degenerate spin- 
ice state and induce various magnetic states— 2 For ex¬ 
ample, with increasing magnetic field along the [100] 
axis, the system enters a Q = 0 state at /Jo H > 0.5 
T— iSl This state has a single spin configuration chosen 
from the degenerate spin-ice ground state, with all the 
spins on each tetrahedron having a component along the 
field direction and actually forming a long-range order. 
When the field is along the [110] direction, the system 
enters a Q = X state at /JoIT > 0.4 T— 2r— The spin 
system is separated into two sets of chains parallel (a) 
and perpendicular (/3) to the field, with long-range fer¬ 
romagnetic order and short-range antiferromagnetic or¬ 
der, respectively— 2r— A more complicated case occurs 
when applying field along the [111] direction. There are 
two successive transitions from the spin-ice state to the 
kagome-ice state and then to the fully-polarized state at 
~ 0.3 and 0.9 T, respectively— L— Another notable char¬ 
acteristics of the spin-ice materials is the very slow spin 
dynamics at low temperatures, demonstrated by the dc 
magnetization, ac susceptibility, and specific-heat mea¬ 
surements, etc i 11 : 14 ! 17 : 29 ^ Due to the slow spin dynam¬ 
ics, these field-induced magnetic transitions show signif¬ 
icant irreversibility, which has been probed by the mag¬ 
netization and neutron scattering ! 14 i 20 i 24 
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Yb 2 Ti 2 07 lias been proposed to be a good candidate 
for the quantum version of spin ice— r— The crystal-field 
structure of the Yb 3+ ion is a ground-state Kramers dou¬ 
blet well-separated from the first excited doublet. The g- 
factor has large planar component, g± = 4.18, compared 
to the [111] component g\\ = 1.77. 37 The net interac¬ 
tion between the neighboring Yb 3+ ions is ferromagnetic 
with a Curie-Weiss temperature of ~ 0.65 K—^ For 
these factors, this material can be described as an effec¬ 
tive pseudospin-1/2 quantum-spin ice, with strong trans¬ 
verse quantum fluctuations of magnetic dipoles—r— The 
nature of ground state of Yb 2 Ti 2 07 is still rather con¬ 
troversial. The specific-heat measurements have re¬ 
vealed a sharp peak at ~ 200 mK, signifying a first- 
order transition. 3 — However, the neutron scattering 
measurements have indicated that the low-T state is 
weakly ferromagnetic, accompanied with short-range and 
dynamic spin correlations. 43 Furthermore, the neutron 
scattering above 200 mK has revealed the pinch-point 
structure, indicating the presence of a spin-liquid phase 
with spin-ice correlations— The excitations of magnetic 
monopoles are expected in this “high”-T quantum spin- 
ice statej 41 ’ 44 ! 45 

Low-temperature heat transport is a powerful tool to 
probe the properties of elementary excitations^— and 
the field-induced magnetic transitions— In principle, 
the magnetic monopoles, as the elementary excitations in 
the spin-ice state, can also contribute to the heat trans¬ 
port by acting as either heat carriers or phonon scatter- 
ers. In a pioneer work by Klemke et al it was con¬ 
cluded that the thermal conductivity (k) of Dy 2 Ti 2 07 
is purely phononic and the field dependence of k is at¬ 
tributed to phonon scattering by monopoles. In contrast, 
Kolland et al. explained their k data of Dy 2 Ti 2 07 in 
terms of the magnetic monopoles making a large contri¬ 
bution to the k in zero field. 6 -~— Toews et al. measured 
the low-T k of another spin-ice material, Ho 2 Ti 2 07 , and 
analyzed the data with the considerations that the mag¬ 
netic monopoles act as heat carriers as well as phonon 
scatterers— Recently, we have also studied the low-T 
thermal conductivity of Dy 2 Ti 2 07 single crystals with 
H || [111] and the heat current parallel or perpendicular 
to the (111) plane. 67 The zero-field k(T) was proved to 
be a purely phononic heat transport, 67 consistent with 
the conclusion by Klemke et al.— Furthermore, an ir¬ 
reversible k(H) behavior was observed in the field re¬ 
gion where the magnetization shows a nearly reversible 
plateau. A picture of the pinning effect of magnetic 
monopoles by the weak disorders was proposed to ex¬ 
plain this novel phenomenon. 6 - It seems that an accurate 
understanding on the low-T heat transport of Dy 2 Ti 2 07 
has not been achieved. 

A very recent work on Yb 2 Ti 2 07 has found a sup¬ 
pression of thermal conductivity in low fields— It was 
explained as a large monopole heat transport at T > 200 
mK. In physics, the quantum monopoles can be more 
promising than the classical monopoles for transporting 
heat because they are dispersive while the classical ones 
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FIG. 1: X-ray diffraction pattern of (110) plane (a and c) and 
the rocking curve of (440) peak (b and d) for two pieces of 
Dy 2 Ti 2 07 and Yb 2 Ti 2 C >7 single crystals, which were orien¬ 
tated by using the X-ray Laue photographs. The full width 
at half maximum (FWHM) of the rocking curves are shown 
in panels (b) and (d). 


are dispersionless. A coherent propagating of monopoles 
was also indicated by the terahertz spectroscopy and 
microwave cavity techniques-® However, more detailed 
studies on the temperature and field dependencies of 
thermal conductivity are desperately needed for estab¬ 
lishing a complete understanding of the heat transport 
properties. 

In this work, we study in details the low-T thermal con¬ 
ductivity of Dy 2 Ti 2 07 and Yb 2 Ti 2 07 single crystals in 
magnetic fields up to 14 T along the [111], [100] and [110] 
directions. It is found that most results of Dy 2 Ti 2 07 are 
not supportive to a sizeable heat transport of monopoles, 
while some results of Yb 2 Ti 2 07 point to a considerable 
quantum monopole transport at T > 200 mK. The k at 
very low temperatures is mainly the phonon transport, 
but the spin fluctuations in Yb 2 Ti 2 07 strongly scatter 
phonons. In both materials, the field-induced magnetic 
transitions can affect the phonon heat transport strongly. 
Moreover, the magnetoelastic coupling are probably also 
playing a role. All these indicate rather complicated 
mechanisms of the heat transport in these materials. 
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II. EXPERIMENTS 

High-quality Dy 2 Ti 2 C >7 and Yb 2 Ti 2 C >7 single crystals 
were grown using a floating-zone technique.' 0 The basic 
low-temperature properties of our crystals, including the 
magnetic susceptibility (down to 2 K) and specific heat 
(down to 0.4 K), were checked and found to be in good 
consistency with most of the reported data from the lit¬ 
eratures. Some of these data were already shown in Ref. 
[roi . To exhibit the crystallinity of samples, some repre¬ 
sentative XRD data on the single crystals are shown in 
Fig. 1. The narrow widths of the rocking curve of the 
(440) Bragg peak (FWHM ss 0.11° and ~ 0.07° for these 
two samples), demonstrates that the crystals have good 
crystallinity. 

For the k measurements, the long-bar shaped sam¬ 
ples were cut from the as-grown crystals along the [111], 
[100] or [110] axes (the uncertainty is smaller than 1°) af¬ 
ter orientation by using back-reflection X-ray Laue pho¬ 
tographs. The k was measured using a “one heater, two 
thermometers” technique in a 3 He refrigerator at 300 mK 
< T < 30 K and a 3 He- 4 He dilution refrigerator at 50 
mK < T < IK, equipped with a 14 T magneti c 60 ! 67 ! 71 
All the measurements were done with both the magnetic 
field and the heat current along the longest dimension 
of the samples, which can minimize the demagnetization 

effects4£ 

The thermal conductivities of Dy 2 Ti 2 C >7 along the 
[111], [100], and [110] axes were measured at low tem¬ 
peratures down to 0.3 K and in magnetic fields up to 14 
T on three samples with the sizes of3.1x0.71x0.15 mm 3 , 
3.08x0.66x0.15 mm 3 and 2.86x0.66x0.15 mm 3 , respec¬ 
tively. In this work, the magnetic-field dependencies of n 
for Dy 2 Ti 2 07 were measured in two different processes. 
The first one is the standard steady-state technique with 
the field changing step by step. 0 ' -60,6 ' That is, the mea¬ 
surements were done by the following steps: (i) change 
field slowly to a particular value and keep it stable; (ii) 
after the sample temperature is stabilized, apply a heat 
power at the free end of sample; (iii) wait some time 
(typically for several minutes at very low temperatures) 
until the temperature gradient on the sample is nearly 
stabilized (judged from the time dependencies of two 
Ru0 2 thermometers on the samples; a relative change 
of less than 1 % within a time window of 2 or 3 min¬ 
utes is the usual criterion); (iv) record the temperatures 
of two thermometers and obtain the temperature gra¬ 
dient. Using this process, the n(H) isotherms for three 
samples were measured at different temperatures and the 
irreversibilities were also probed with the measurements 
done in changing field (from zero after zero-field cooling) 
up and down in the step mode. The speed of changing 
field is slow enough to avoid observable heating effect 
from the eddy current. More exactly, the field ramp¬ 
ing rate is not a constant, and it must be very slow at 
HoH < 0.1 T (0.001-0.01 T/min) and could be quicker 
at HoH > 0.1 T (0.01-0.05 T/min). Note that this pro¬ 
cess is suitable for most of solid materials. However, it is 


known that Dy 2 Ti 2 07 has very slow spin dynamics at low 
temperatures^ which may cause some relaxation phe¬ 
nomenon in the temperature gradient measurements ) 62 ’ 64 
if the magnetic excitations are involved in the heat trans¬ 
port behaviors. For this reason, the n measurement done 
with the sweeping-field mode could not be very mean¬ 
ingful because the thermal relaxation and spin relaxation 
must be strongly mixed. 

To probe the time relaxation of the k, the measure¬ 
ments were also done in a second process, which is dif¬ 
ferent from the first one in the last step. Namely, after 
the field and sample temperature are stabilized, apply a 
suitable heat power at the free end of sample and wait 
for about 30 minutes; then, record the time dependence 
of temperature gradient on the sample. Since the relax¬ 
ation was found to be extremely slow at some particular 
fields, it is not possible to finish a relaxation measure¬ 
ment at different fields in a up and down loop. The relax¬ 
ations were measured for three different samples at very 
low temperatures. In this case, the k at a certain field 
was measured separately in the field-increasing and field- 
decreasing processes, respectively. In the former case, the 
sample was cooled to the measurement temperature (for 
example, 0.36 K) in zero field and then field was changed 
slowly to the particular values, and then the k was mea¬ 
sured in each fields. In the latter case, the sample was 
first cooled to 0.36 K in zero field; after that the field 
was swept slowly to 1-2 T, which are well above the irre¬ 
versible region as we show below, then decreased slowly 
to the target field. 

The thermal conductivities of Yb 2 Ti 2 07 along the 
[100], [110], and [111] axes were measured at low tem¬ 
peratures down to 50 mK and in magnetic fields up to 14 
T on three samples with dimensions of2.07x0.64x0.12 
mm 3 , 2.87 x 0.57 x 0.14 mm 3 and 2.53 x 0.60 x 0.14 
mm 3 , respectively. The measurements of Yb 2 Ti 2 07 sin¬ 
gle crystals are much simpler, since there is no observable 
relaxation effect. 


III. RESULTS 

A. Zero-Held k(T) of Dy^THCU 

Figure 2 shows the temperature dependencies of n for 
Dy 2 Ti 2 07 along three directions and in zero field. Appar¬ 
ently, the heat transport is nearly isotropic and shows a 
simple phonon transport. The phonon peaks are located 
at ~ 15 K and a rough T 2 " 5 dependence of n is visi¬ 
ble at subkelvin temperatures, which is however weaker 
than the standard T 3 behavior of phonon thermal con¬ 
ductivity at the boundary scattering limit 4^ Similar re¬ 
sult has been obtained in the earlier works and the mag¬ 
netic scattering on phonons was discussed to be impor¬ 
tant at temperatures below 10 K.—x 65 - 67 The mean free 
path of phonons at low temperatures are calculated in 
the usual wayj 58,6 ' - ' 2 The phononic thermal conductivity 
can be expressed by the kinetic formula K p h = Cv p l/3r - 
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FIG. 2: (Color online) Temperature dependencies of the ther¬ 
mal conductivities of Dy 2 Ti 2 07 in zero field for heat current 
along the [111], [100] or [110] axis. The dashed line shows the 
T 2 ' 5 temperature dependence. The inset shows the tempera¬ 
ture dependencies of the phonon mean free path l divided by 
the averaged sample width W for these samples. 

where C — /3T 3 is the phonon specific heat at low tem¬ 
peratures, v p is the average velocity and l is the mean 
free path of phonons. The j3 value of Dy 2 Ti 2 07 crystals, 
obtained from our specific-heat data,. 70 ' is 1.14 x 10~ 3 
J/K 4 mol (which is in good consistency with the result 
from other group2£). The inset to Fig. 2 shows the ratios 
l/W for three samples, where W is the averaged sample 
widthi 58 ! 67 ' 72 It is clear that the ratios increase quickly at 
very low temperatures and are expected to approach 1 at 
T < 0.3 K. This calculated result is compatible with the 
temperature dependence of n, which is a bit weaker than 
T 3 . Note that if there were other type of heat carriers 
(e.g. magnetic monopoles) that makes a large contribu¬ 
tion to the k. the actual phonon mean free path would be 
much smaller than those in Fig. 2, which is apparently 
not very reasonable. 


B. k(H) of Dy2Ti 2 C>7 

The low-T thermal conductivities of Dy 2 Ti 2 0y had 
been studied in magnetic fields up to only 7 T and the 
main finding was that the k can be strongly suppressed 
by rather weak fields (< 1 T)^ ~ 65 i 67 There are two dif¬ 
ferent explanations for the field-induced suppression of n 
in Dy 2 Ti 2 C> 7 . Namely, it is due to either the suppression 
of a monopole heat transport or the field-enhanced mag¬ 
netic scattering on phonons£ 2 ~ 65 ’ 6 £ Note that the k(T) 
data shown in Fig. 2 are not supportive for a sizeable 
monopole heat transport in zero field. In passing, a re¬ 



FIG. 3: (Color online) Magnetic-field dependencies of thermal 
conductivities of Dy 2 Ti 2 07 : (a,c,e) data in fields up to 14 T 
and along three directions; (b,d,f) the zoom-in of low-field 
data. In this figure, the magnetic field is increased step by 
step after cooling the samples in zero field. The directions of 
the magnetic field and heat current ( Jh) are the same. 


cent work on another spin-ice material, Ho 2 Ti 2 07 , pro¬ 
posed that magnetic monopoles are playing a dual role in 
the heat transport; that is, they can both transport heat 
and scatter phonons^ In that work, the contribution of 
the monopole transport to the total thermal conductiv¬ 
ity was estimated to be about 10 %*££ much smaller than 
that proposed for Dy 2 Ti 2 07 ^— 

In present work, the magnetic-field dependencies of k 
are measured up to 14 T, as shown in Fig. 3. The low- 
field behaviors are almost the same as those results in 
earlier studies. For H || [111], at very low temperatures, 
the k display two anomalies in low fields: at 0.36 K, the 
k(H) curve shows a quick decrease and a peak-like tran¬ 
sition at about 0.25 and 0.9 T, respectively, which are in 
good correspondence with the critical fields of the sub- 
sequential transitions from the low-field spin-ice state to 
the kagome-ice state then to the saturated state (“3-in, 
1-out” or “3-out, 1-in”) ^22 At 0.5 K, the k(H) curve 
behaves like two step decreases at these two transitions. 
The sharpness of the transitions is quickly smeared out 
upon increasing temperature. For H || [100], the 0.36- 
K k{H) curve shows a steep and broad step-like transi- 
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tions at 0-0.2 T and 1-2 T, respectively. Upon increas¬ 
ing temperature, the first transition gradually becomes 
weaker while the second one is enlarged and shifts to 
higher fields. At 5 K, the first step is completely gone 
and the second one evolves into a broad transition at 1-4 
T. The first transition is in good correspondence with the 
critical field of a long-range-order transition^ in which 
the external magnetic field chooses one of the degenerate 
ground states with all the spins having the components 
along the field. For H || [110], the low-field k(H) behave 
rather similarly to the case of H || [100]. For example, at 
0.36 K the k(H) curve also shows two step-like decreases 
at 0-0.25 and 0.8-1.5 T. The first step-like change of k 
at 0.25 T is in good correspondence with a first-order 
transition from the low-field spin-ice state to the spin- 
chain state^2r— in which the a-chains form a long-range 
ferromagnetic order and the ^-chains form a short-range 
antiferromagnetic order. Upon increasing temperature, 
the first transition gradually becomes weaker, while the 
second one becomes larger and shifts to higher fields. 

The drastic change of the heat transport at the mag¬ 
netic phase transition has been widely observed in many 
magnetic materials and is usually related to the evolu¬ 
tion of magnetic excitations—However, the particu¬ 
lar features of k(H) at the magnetic transitions can be 
significantly different from each other, which depend on 
both the natures of magnetic transitions and the roles 
of magnetic excitations in the heat transport (as heat 
carriers or phonon scatterers). In the spin-ice state of 
E^TUCU, the monopoles are the only one type of mag¬ 
netic excitations. In this regard, it was naively assumed 
that the quick decreases of k at the low-field magnetic 
transitions are caused by the suppression of monopole 
transport—This scenario will be further analyzed in 
the following sections. 

A remarkable finding is that these n{H) curves still 
exhibit clear or strong field dependencies when the field 
is above 7 T, which have not been achieved in those ear¬ 
lier works—~ 65,6 ' At subkelvin temperatures, the n(H) 
for H || [111] and [100] show a broad-valley-like feature; 
that is, after arriving a minimum at several tesla, the n 
gradually increases in high magnetic fields. In particular, 
at 0.36 K, the k at 14 T || [100] is already larger than that 
at 0.5 T, where the thermal conductivity was taken for 
the phononic background in a recent workj 63 In addition, 
the k(H) for H || [111] and [100] are not yet saturated in 
14 T. 

It is notable that the high-field n(H) behave rather 
differently for H || [110]. The k(H) show weak field de¬ 
pendence at several tesla, particularly at 0.36 K, which 
is also consistent with the earlier reported data in field 
up to 7 T— However, at low temperatures the k show 
another strong decrease at fioH > 7 T and there is no 
sign that the n can be recovered even at 14 T. Since this 
phenomenon does not appear for fields along other di¬ 
rections, a high-field polarization effect is not likely the 
reason. For now, there are no other experimental results 
for Dy 2 Ti 2 C >7 in [110] field higher than 7 T, so it is un¬ 


known whether the strong suppression of n at fJ-oH > 7 
T (|| [110]) is related to some magnetic transitions. For 
clarifying, high-field investigations by using other mea¬ 
surements are necessary in the near future. 

The paramagnetic scattering effect related to the 
crystal-field levels of Dy 3+ ions was proposed to play an 
important role in the field dependencies of k i 61 i 67 It is 
known that the Dy 3+ ions have a degenerate doublet of 
the lowest crystal-field level. In magnetic field, a Zeeman 
splitting of the doublet can produce resonant scattering 
on phonons and give both a low-field suppression and a 
high-field recovery of k, as some other magnetic materi¬ 
als have shown—— — However, since the paramagnetic 
scattering is a qualitatively isotropic effect, the obvious 
anisotropic and complicated high-field behaviors in Fig. 
3 should not be solely relevant to this mechanism. 

C. Irreversibilities of k(H) of Dy2Ti2C>7 

It is found that the magnetic-field dependencies of k 
display irreversible behavior at very low temperatures for 
all three field directions, as shown in Fig. 4. These data 
reproduce our earlier results for H || [111] and are similar 
to those reported by other group^ 2-65,6 ' The reversibil¬ 
ity weakens quickly with increasing temperature and are 
hardly visible at T > 0.5 K. Note that this is a bit dif¬ 
ferent from some earlier results that clear hysteresis was 
observed at 0.6 K—— The reason is that those data 
were taken in a sweeping-field mode. 

There are significant differences in the irreversibility 
between the k(H) data with H || [111] and the data taken 
with H || [100] and [110]. First, the hysteresis of k{H) 
with H || [111] is located mainly in a field region be¬ 
tween the two field-induced magnetic transitions, while 
those with H || [100] and [110] are located at the magnetic 
transitions. Actually, the irreversibility for H || [100] and 

[110] seems to be a common “supercooling” phenomenon 
of a first-order phase transition. Note that the ultra-low- 
T neutron scattering also demonstrated clear hysteresis 
of the field dependencies of Bragg scattering intensity for 
H || [100] and [110]— Second, the irreversibility for H || 

[111] appears in the kagome-ice state, in which the mag¬ 
netization curves do not show any hysteresis ; 13 i 24 i 62 i 63 i 67 
whereas those for H || [100] and [110] have good corre¬ 
spondence with the magnetization hysteresis — — Third, 
the magnitudes of k are larger for decreasing fields than 
those for increasing field with H || [111], while the oppo¬ 
site behavior is observed with H || [100] and [110]. All 
these clearly indicate that the origin of the irreversibility 
with H || [111] is different from others. 

D. Relaxation phenomenon of k of Dy2Ti2C>7 

It is known that the thermal conductivity of Dy 2 Ti 2 C >7 
shows not only irreversible behaviors but also a time 
relaxation—-® 4,67 In present work, the relaxation effect 
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FIG. 4: (Color online) The low-T k(H) loops with magnetic 
field sweeping up and down along the [111], [100], or [110] di¬ 
rections. The data shown with solid symbols are measured in 
ascending field after the sample is cooled in zero field, while 
the open symbols show the data with descending field, as 
also indicated by the arrows. These data were measured in 
the normal way, that is, the waiting time for the tempera¬ 
ture stabilization is in a scale of ten minutes. Note that the 
irreversibilities are not visible at T > 0.5 K (not shown here). 


is carefully studied and the representative data are shown 
in Fig. 5. Note that at 0.36 K, the relaxation of k can be 
very strong and very slow, which prevents us from get¬ 
ting the equilibrium state in the available time window 
of our 3 He refrigerator. However, the relaxation effect is 
not visible in high magnetic fields or at T > 0.5 K, where 
the k(H) hysteresis in Fig. 4 disappears. 

Among three field directions, the relaxation of k in H || 
[111] is the weakest. As shown in Figs. 5(a) and 5(b), in 
either zero held or 0.5 T (the Held showing largest hys¬ 
teresis), the magnitude of relaxation is very small. Typ¬ 
ically, after waiting for about 30 min (the starting point 
of the recorded data in Fig. 5), the measured k change 
in a scale of only 1 -2 %, which is nearly the same as 
the error of this measurement. This result confirms our 
earlier work that the relaxation effect in H || [111] is not 
significant. 6 ' Although one cannot completely rule out 
the possibility of a relaxation with very long time scale 
(e.g. much longer than the time window of this mea- 
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FIG. 5: (Color online) Time dependencies of thermal con¬ 
ductivities for field along the [111], [100], or [110] directions 
and at 0.36 K. The details for these measurements are de¬ 
scribed in the experimental section. The open squares show 
the data measured in ascending field after the sample is cooled 
in zero field, and the open circles with descending field from 
high field, where no relaxation and irreversibility are observed. 
The horizontal lines are guides for eyes while the curves are 
fittings using formula ©• All the k( 0) are taken the values 
from the normal measurements, as shown in Figs. 3 and 4. 


surement), the magnitude of the k changing with time 
for if || [111] is significantly smaller than those for other 
Held directions. A bit larger relaxation of k for H || [111] 
could be observed in a different measurement process, in 
which the magnetic field was continuously swept, 64 
The relaxation of k is the most significant for H || [100]. 
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As shown in Fig. 5(c), although there is no sizeable re¬ 
laxation effect after cooling in zero field, decreasing field 
to zero from high field (1 T or higher) results in a slow 
relaxation. This is very different from the case of [111] 
field shown in Fig. 5(a). In addition, in non-zero fields, 
the relaxation could be very strong and is related to the 
history of applying field. For example, increasing field 
from 0 to 0.1 T does not cause significant relaxation, but 
decreasing field from high value to 0.1 T leads to a very 
large relaxation; the n can be changed about 10 % and is 
not saturated after a very long time of about 2400 min, 
as shown in Fig. 5(d). A contrary feature is observed at 
0.15 T; as shown in Fig. 5(e), decreasing field from high 
value to 0.15 T does not cause significant relaxation, but 
increasing field from 0 to 0.15 T leads to a very slow relax¬ 
ation. Furthermore, the k slowly relaxes towards larger 
values in the case of decreasing field to 0.1 T, while it 
slowly relaxes towards smaller values if increasing field 
to 0.15 T. Such differences demonstrate that the hystere¬ 
sis of k(H) is directly related to the relaxation effect, as 
discussed in the following text. 

The time dependence of k can be fitted by a single 
exponential function, 


«(#) 

«(o) 


= a 


be~ t/T , 


( 1 ) 


where r is the characteristic time constant, and a and b 
are time independent parameters. The time constant r 
could be as large as 1100-1600 min. 

The situation for H || [110] is very similar to that for 
H || [100], but the relaxation is a bit weaker. For ex¬ 
ample, the time constants r are about 900 min and the 
k changes about 5 % with decreasing held to 0.1 T, as 
shown in Fig. 5(f). 

The slow spin dynamics of the spin ice has been 
probed by other experiments on Dy 2 Ti 2 C> 7 , such as 
the dc magnetization, ac susceptibility, and specific 
heat i 11 1 14 i 17 i 29 ~ — It has been found that the spin relax¬ 
ation time in the spin-ice state showed an exponential 
increase at low temperatures and is greater than 10 4 s 
at T < 0.45 K. The current explanations are mainly in 
terms of diffusive motion of monopoles. After cooling 
in zero held to a low but hnite temperature, the density 
of monopoles decreases gradually with time. This could 
induce a slow relaxation with large time scale. How¬ 
ever, it is notable in the present k measurements that 
cooling in zero held does not lead to a signihcant re¬ 
laxation of k. This means that the monopole dynamics 
does not have close correlation with the heat transport of 
Dy 2 Ti 2 C> 7 . On the other hand, the magnetization mea¬ 
surements with the held-quench or pulse-held processes 
have revealed the monopole dynamics having a time scale 
of only several minutes This is in contrast with the 
extremely slow relaxation of n in magnetic held. More 
importantly, if the relaxation of n in hnite helds, as shown 
in Figs. 5(d- f), were mainly caused by the monopole dy¬ 
namics (either transporting heat or scattering phonons), 
one would expect opposite trends of n(t) for increasing 



FIG. 6: (Color online) The 0.36-K k(H) loops with the mag¬ 
netic field sweeping up and down along the [111], [100], or 
[110] directions. The data shown with squares are measured 
in the normal way, that is, the waiting time for the tempera¬ 
ture stabilization is in a scale of several or ten minutes. The 
circles show the data extrapolated at t —» oo from the time- 
relaxation measurements. The arrows indicate the direction 
of changing held. 


or decreasing to the same held. Again, none of the data 
supports this expectation. Therefore, the relaxation of k 
cannot be well explained based on the assumption that 
the monopoles play an important role in the heat trans¬ 
port properties. Furthermore, the anisotropic relaxation 
of k for different held directions indicates that the slow 
spin dynamics associated with the held-induced transi¬ 
tions for H || [100] and [110] are strongly coupled with 
phonons, whereas that for H || [111] is not. 

The relaxation httings using formula m can get not 
only the time constant but also the extrapolated values 
of k(H)/k( 0) at t —> oo. Then, the magnetic-field depen- 











dencies of n in the expected “equilibrium or completely 
relaxed state” can be obtained. Figure 6 shows the com¬ 
parison of the k(H) between the data of normal measure¬ 
ments and the data after complete relaxation for all three 
field directions. The main result is that the irreversibil¬ 
ity in the [111] field changes only slightly and still clearly 
exists after long-time relaxation, but the irreversibilities 
in the [100] and [110] fields almost disappear for the data 
extrapolated at t —> oo. This demonstrates that the ir¬ 
reversibilities of k{H) in the [100] and [110] fields are re¬ 
lated to the extremely slow spin dynamics in Dy 2 Ti 2 C> 7 , 
while that in the [111] field has a different origin. As dis¬ 
cussed in our previous work, a pinning effect of magnetic 
monopoles by crystal disorders can explain the particular 
irreversibility with H || [111]. 67 It is also notable that in 
this case the irreversibility of k(H) is actually not large, 
which is consistent with the previous conclusion on the 
minor role of monopoles in the heat transport. 


E. Analysis on the Dy2Ti2C>7 data assuming the 
monopole transport 

The present experimental results have already met 
with some difficulties in explaining the data when as¬ 
suming large monopole heat transport. A quantitative 
analysis on the low-field decrease of n is useful to make 
further clarification. Figure 7 shows the schematic pic¬ 
ture of monopole motions that can contribute to the heat 
transport along the magnetic-field directions. 

The changes of energy gap between the ground state 
and excited state and the density of monopoles are depen¬ 
dent on the held direction. 68 ' 76 When the held is applied 
along the [100] direction, the ground state and excited 
state split into three and two levels, respectively^ 6 Be¬ 
cause the split levels are rather complicated, a simplifica¬ 
tion should be carried out. According to the calculation 
of Castelnovo et ol.r^ the density of monopoles decays in 
an approximately exponential speed because the energy 
gap increases. Here, only the hopping from the lowest 
level of ground state to the lowest level of excited state 
is considered, whose gap is A = A 0 + (2/v / 3 )g^sH and 
Ao = 4.35 K is the energy gap in zero held. 76,77 Figure 
7(a) presents this process, in which the spin 1 hips and 
a pair of monopoles m a and mj is created, mb has three 
choices of spin hip to recover to the ground state. If the 
spin 1 hips, the two monopoles annihilate, and there is 
no monopole movement. If one of the other two equiva¬ 
lent spins 2 or 4 hips, mb moves to the next tetrahedron 
and becomes m c , as shown by Fig. 7(b). This move¬ 
ment has a positive component along the heat current 
and makes contribution to thermal conductivity. After 
this process involving the heat transport, the energy of 
the two tetrahedra is lifted by E = (2/\/3On 
one hand, increasing held leads to larger E. As a re¬ 
sult, this process becomes more and more difficult and 
the movement of monopoles are restricted. On the other 
hand, the energy gap A for monopoles also increases and 



FIG. 7: (Color online) Schematic plot for the motions of mag¬ 
netic monopoles under magnetic fields along three directions. 
Note that only those motions that can contribute to heat 
transport in the [100] (a,b), [110] (c,d) and [111] (e-h) direc¬ 
tions ( Jh || H) are shown. Red and blue dots represent mag¬ 
netic monopoles with opposite polarities. Blue arrows repre¬ 
sent the flipped spins, associated with a monopole movement 
from the dashed-dot position to the solid-dot position. 


fewer monopoles can be excited. 

In helds along the [110] direction, both the ground 
state and the excited state split into three levels. With 
the same simplification, the energy gap between the low¬ 
est levels of ground state and excited state is A = A 0 ; 
that is, the gap for monopoles is field independent for H || 
[110]. This means the density of monopoles also does not 
vary with field. Figure 7(c) shows a pair of monopoles, 
m a and mb, created by this process. Apart from the an¬ 
nihilation caused by the spin 3, the inequivalent spins 2 
and 4 provide other two recovery choices for mb- If the 
spin 4 flips, the monopole propagates along the direction 
perpendicular to heat current, and this process makes 
no contribution to heat transport. If the spin 2 flips, 
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as shown by Fig. 7(d), mb moves to the next tetrahe¬ 
dron and becomes m c , with a positive component along 
the heat current. This process accompanied by the heat 
transport leads to an energy lift of E = (2-\/6 /3)gpsH. 
This energy lift increases with field, and therefore the 
monopole movement is suppressed. 

In the case of if || [111], the propagation of monopoles 
along the [111] direction is the most complicated. The 
ground state and excited state split into two and four 
levels in the [111] field, respectively. The excitation gap 
of monopoles m a and mb, as shown in Fig. 7(e), is 
A = Ao — (2/3 )gpgH, which results in an increase of 
monopole density. For the movement of mb, two effec¬ 
tive flips caused by the spins 2 and 4 are equivalent and 
lift the energy by E\ = (2/3)g/isif. However, m c , which 
is generated in this process shown by Fig. 7(f), does not 
leave the original (111) plane. Further effective flips have 
to be considered. First, if the next flipping occurs with 
the spin 3, as shown in Fig. 7(g), m c moves within the 
(111) plane and becomes nid . Though energy drops by 
E 2 = (2/3 )gpsH in this process, no heat transport along 
the [111] direction occurs. Second, if the spin 1 flips, it 
can be seen in Fig. 7(h) that m c goes across the (111) 
plane and contributes to heat transport. m c turns into 
m e and energy is lifted by £3 = 2 gpBH. Thus, two- 
step flips are needed to transporting heat, and the total 
energy lift £ = £1 + £3 also increases with field. 

The above discussions indicate that magnetic field 
affects both the excitation and the propagation of 
monopoles. The former refers to the monopole energy 
A and the density p , while the latter refers to veloc¬ 
ity v. As discussed above, the monopole energy A with 
H || [100], [110] and [111] are A 0 + (2/v / 3 )gpsH, A 0 and 
A 0 — (2/3 )gp,BH, respectively, as shown in Fig. 8 (a). For 
the monopole density, only free monopoles are concerned. 
According to the Debye-Hiickle theory^- the monopole 
density is determined by the energy gap, 

2 exp (A/fcflT) 

1 + 2 exp (A/ksT) 

This formula gives the field dependence of the monopole 
density, as shown in Fig. 8 (b). The specific heat of 
monopoles can be obtained by 2 i 

Cm = N-^(pA), (3) 

where N represents the number of tetrahedra per unit 
volume. Since the positive and negative monopoles 
would annihilate when they meet, the mean free path 
of monopoles is assumed to be the averaged distance of 
monopoles with the same polarity. Thus, it can be writ¬ 
ten as 


lm = {N P -)-C\ ( 4 ) 

If the velocity of monopoles is independent of field, 
according to the kinetic formula i^rn — m^m/3, 



FIG. 8: (Color online) (a) Magnetic-field dependencies of 
the energy gap of monopole excitations in Dy 2 Ti 207 . (b) 

Magnetic-field dependencies of the monopole density, (c) Cal¬ 
culated monopole thermal conductivity using the kinetic for¬ 
mula after considering the field dependencies of A, p and ve¬ 
locity of monopole. The inset shows the calculation using the 
same formula, but the field dependencies of velocity are not 
considered, (d) The monopole thermal conductivity extracted 
from the experimental data in Fig. 6, assuming the decrease 
of raw-data k from 0 to 0.5 T is just a decrease of K m . All 
the calculations and experimental data are taken at 0.36 K. 


K m {H)/K m (0) at 0.36 K can be calculated and are shown 
in the inset to Fig. 8 (c). This calculation shows that the 
K m increases strongly with H || [111] and is independent 
of field with H || [110], which are completely different 
from the experimental observations. In fact, the velocity 
of monopole should be sensitive to the magnetic field and 
its direction, as discussed above. Currently, the velocity 
of monopoles can not be calculated accurately. If only 
the relative change is concerned, based on the additional 
energy lift E caused by the magnetic field, we simply 
assume that the monopole velocity decays exponentially 
with increasing E; that is, 


Vm 


v m0 exp ( 


E 


aksT 


(5) 


where v m o is the monopole velocity at zero field and 
a is an adjustable parameter. E are given by £ = 
( 2 /V3)gp B H, E = (2v / 6/3 )gp B H and E = (8/3 )gp B H 
for H || [100], [110] and [111], respectively. With this 
additional assumption, n m (H)/K m ( 0) at 0.36 K are cal¬ 
culated and shown in Fig. 8 (c). 

For comparison, Fig. 8 (d) shows the low-field “« m 
data” at 0.36 K extracted from the experimental results 
for H || [100], [110] and [111]. Here, the “relaxed” k(H) 
data from Fig. 6 are used, since the above calculations do 
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not take into account the relaxation effect. It is assumed 
that the decrease of k from 0 to 0.5 T is simply due to 
the suppression of K m . The raw-data of k at 0.5 T is 
taken as the background from the phonon transport. It 
seems that there is a rather good consistency between 
the calculations and the experimental data. 

Furthermore, from the assumed K m data shown in Fig. 
8 (d), the absolute values of K m (0) along the [100], [110] 
and [111] directions are 2.39 x 10 -2 , 8.85 x 10 -3 and 
6.01 x 10 ~ 3 W/Km, respectively. This anisotropy can 
be explained as the difference in monopole ratio, n, that 
conducts heat in different directions. Based on the above 
discussion, only certain spin flips contribute to transport¬ 
ing heat. If we assume the probabilities of spin flips in 
one tetrahedron are equal and omit flips leading to all-in 
or all-out state, the probability for a monopole conduct¬ 
ing heat is 2/3, 1/3 and 2/9 in the [100], [110] and [111] 
directions, respectively, which gives the monopole ratio, 
n, participating in transporting heat. From Eqs. (2), (3) 
and (4), the monopole density, specific heat and mean 
free path at 0.36 K are 1.13 x 10 -5 per tetrahedron, 
176.88 J/Km 3 and 2.83 x 10 -8 m, respectively. Using 
the formula K m = nC m v m ol m /3, the v m o in the [100], 
[110] and [111] directions can be calculated to be 2.12 x 
10 4 , 1.59 x 10 4 and 1.62 x 10 4 m/s, respectively. These 
values are unreasonably too large since monopoles are 
known to be dispersionlessi!! 

Finally, we need to point out another experimental 
phenomenon that cannot be explained using the assump¬ 
tion of monopole heat transport. As found in our former 
work, for H || [111] the field dependencies of k along 
the field and perpendicular to it are nearly isotropic**! 
As discussed above, these two perpendicular monopole 
transport are established with the processes shown by 
Figs. 7(e)-7(f) 7(h) and Figs. 7(e)-7(f)-7(g), respec¬ 
tively. In the former case, the process in Fig. 7(h) is 
suppressed most strongly by magnetic field. According 
to magnetic susceptibility result, the spin 1 in Fig. 7(e) 
is fixed along held direction at about 0.3 Twhich leads 
to the disappearance of K m (H). In the latter case, mag¬ 
netic held only suppresses the process in Fig. 7(f). At 0.3 
T, although hxing spin 1 does not forbid the monopole 
transport directly, it reduces the degree of freedom of 
monopoles and forbids the subsequent monopole propa¬ 
gating after the hip of spin 1. Therefore, n m {H) should 
have a kink-like decrease at 0.3 T instead of completely 
disappears. When held reaches 0.9 T, the spin 2 in Fig. 
7(f) is also fixed!! and the K m is completely suppressed. 
Therefore, the K m (H ) along the [111] held and perpen¬ 
dicular to it are expected to present obvious anisotropy, 
which is not supported by the experimental results**! 

The analysis and calculations indicate that only some 
of the low-held results of k(H) can be qualitatively ex¬ 
plained as the suppression of monopole thermal con¬ 
ductivity in magnetic fields. However, the quantita¬ 
tive analysis also points out that the monopoles are not 
likely making a large contribution to the heat trans¬ 
port. It is essentially coincided with some other exper- 



T (K) 


FIG. 9: (Color online) Temperature dependencies of the ther¬ 
mal conductivity along the [ 100 ], [ 110 ] and [ 111 ] axes of 
Yb 2 Ti 2 C >7 single crystals in zero field. The dashed line shows 
the T 2 dependence which the k follow between 200 mK and 
2 K. The inset shows the temperature dependencies of the 
ratio of the phonon mean free path l to the averaged sample 
width W. The l is calculated assuming that the k is purely 
phononic. 


iments. For example, it has been found that the hop 
rate of monopoles is roughly proportional to the density 
of monopoles at low temperatures (< 1.5 K)j!£ Further¬ 
more, a recent adiabatic susceptibility measurement re¬ 
vealed a Brownian characteristic of the monopole move¬ 
ment in the spin-ice state*!! with a much smaller diffu¬ 
sion constant than that proposed by some heat transport 
work . 62 There should be more important factors that are 
responsible for the field dependencies of k in Dy 2 Ti 2 07 . 


F. Zero-field k(T) of Yb 2 Ti 2 07 

Figure 9 shows the temperature dependence of n of 
Yb 2 Ti 2 C >7 single crystals in zero field. Similar to the 
case of Dy 2 Ti 2 07 , the thermal conductivity is nearly 
isotropy for the heat current along the [ 100 ], [ 110 ] and 
[111] axes. The peak at 13 K is the typical feature of 
a phonon transport in insulators. Moreover, there are 
two notable features of these data. First, the tempera¬ 
ture dependencies display kink-like anomalies at ~ 200 
mK where the slopes of k(T ) changes strongly. This 
is likely related to the first-order transition of mag¬ 
netism. Second, at subkelvin temperatures, the k do not 
show the T 3 dependence, which is expected for phonon 
thermal conductivity in the boundary scattering limit. 
Actually, the n nearly follow the T 2 dependence from 
about 2 K to 200 mK. This is consistent with a recent 
report;!! Although the well-known T 3 ballistic behavior 
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of phonons has been rarely observed in the transition- 
metal compounds, including the high-T c cuprates, the 
multiferroic manganites and the low-dimensional quan¬ 
tum magnetS)^ ~ 51 ’ 53 i 55 i 57 i 58 ’ 78 i 79 the k can usually ex¬ 
hibit the temperature dependence rather close to T 3 . 
In Dy 2 Ti 2 C> 7 , as shown in Fig. 2, the T 2 ' 5 dependence 
of k is observed at temperatures of several hundreds of 
millikelvins. In Yb 2 Ti 2 C> 7 , however, the T 2 dependence 
differs too much from the T 3 behavior. There are two 
possible reasons for this deviation from the boundary 
scattering limit. If there are sizeable scattering effects 
on phonons even at the lowest temperature, the phonon 
mean free path is significantly smaller that the sample 
size and T 3 behavior cannot be observed. Another pos¬ 
sibility is the presence of the other term of heat conduc¬ 
tivity, for example, from the magnetic excitations. In 
any case, the estimation of phonon mean free path can 
provide some useful information. 

With the /? value (= 1.61 xlO~ 3 J/K 4 mol) from the 
specific-heat data^ 9 the phonon mean free path can be 
calculated assuming the k is purely phononic. The in¬ 
set to Fig. 9 shows the temperature dependencies of 
the ratio l/W. It is found that the l/W ratios are only 
about 0.3-0.4 at the lowest temperatures, indicating that 
the boundary scattering limit of phonons is not estab¬ 
lished and there are still some microscopic scattering of 
phonons. It should be noted that the above calculation is 
based on the assumption of a purely phononic heat trans¬ 
port. If there were a magnetic term of heat transport, the 
phonon mean free path would be even smaller. Neverthe¬ 
less, it is clear that the phonon scattering is stronger in 
Yb 2 Ti 2 07 , compared to Dy 2 Ti 2 07 . It is understandable 
from the magnetic scattering since Yb 2 Ti 2 07 is known 
to have strong spin fluctuations^ 

All the zero-field k(T) curves show a kink-like decrease 
at 200 mK with lowering temperature. It might be that 
some term of thermal conductivity disappears when en¬ 
tering the low-T state. This is likely a signature that the 
magnetic excitations (monopoles) in the quantum spin- 
ice state (at T > 200 mK) can act as heat carriers: 68 


G. k(H) of Yb 2 Ti 2 0 7 

Figure 10 shows the magnetic-field dependencies of k 
at low temperatures for these Yb 2 Ti 2 07 single crystals. 
The data have two remarkable features. First, with in¬ 
creasing field, the k exhibits strong changes in low fields, 
followed by a quick increase, and finally becomes field in¬ 
dependent in high fields. It naively indicates that there 
are magnetic scattering effects on phonons that can be 
suppressed in high fields, where the spins are polarized. 
This is consistent with the zero-field k{T) data and the 
small phonon mean free path. The strong spin fluctu¬ 
ations revealed by the neutron scattering is most likely 
the source of phonon scattering in zero field. Second, 
the low-field behaviors are somewhat different for differ¬ 
ent directions of magnetic field. For H || [100], the k 




FIG. 10: (Color online) Magnetic-field dependencies of the 
thermal conductivity of Yb 2 Ti 2 0 7 crystals at low tempera¬ 
tures. Panels (a,c,e) show the data in fields up to 14 T and 
along three directions, while panels (b,d,f) show the zoom-in 
data of the low-field plots. Since the magnetic fields are ap¬ 
plied along the direction of the heat current (Jh), the demag¬ 
netization effect is negligibly small for these long-bar shaped 
samples. 


increase monotonically with field when T < 200 mK, but 
show a low-field valley when T > 200 mK. The position 
of dip field shifts to higher field with increasing temper¬ 
ature. For H || [110] or [111], the k{H) isotherms at 
the lowest temperatures show a “dip” at 0.5 T. In par¬ 
ticular, the dip is very sharp for H || [110], reminiscent 
of a field-induced transition. In this regard, the absence 
the “dip” for H || [100] indicates that either there is no 
field-induced transition or the transition field is very low 
(smaller than 0.1 T, which is the step size of our n(H) 
measurements). 

As suggested in a recent work, the low-field decrease 
of k at T > 200 mK may be due to the suppression of 
monopolc heat transport in magnetic field. 68 However, it 
should be pointed out that the mechanism of monopole 
transport is far from well understood. The present k(H) 
data with field along three different directions show some 
peculiar phenomenon. To illustrate, at T > 200 mK 
the low-field behavior of k are qualitatively different for 
different field directions. The position of dip field shifts 
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to higher field with increasing temperature for H || [100], 
while it changes very weakly for H || [110] and [111]. This 
anisotropic behavior calls for a more precise theoretical 
description of the monopole transport. 

Qualitatively, a low-field valley and a strong enhance¬ 
ment at high field in k(H) are also similar to the para¬ 
magnetic scattering effect i’ 59 ’ 73- — It is known that the 
Yb 3+ ions have a Kramers doublet of the lowest crystal- 
field level. In magnetic fields, a Zeeman splitting of the 
doublet can produce resonant scattering on phonons and 
give both a low-field suppression and a high-field recov¬ 
ery of k, 59 i 73 ~ — A simplified calculation has shown the 
qualitative behaviors of n(H) in the case of paramagnetic 
scattering; one important feature is that the position of 
the k(H) minimum shifts to higher field with increas¬ 
ing temperature. 74 ' 75 Since the changes of low-field valley 
with increasing temperature are not in good consistency 
with this calculation, the paramagnetic scattering may 
not be important in Yb 2 Ti 20 7 . 

The low-field behavior of k at T < 200 mK displays dif¬ 
ferent characteristics from the data at T > 200 mK. The 
low-T magnetization measurements have indicated that 
the ground state is ferromagnetic with the spin easy-axis 
in the [100] direction, but the residual spin fluctuations 
are evidenced^ These results indicated that the spins 
are easily polarized in a rather weak field. In this regard, 
the quick increase of k with field seems closely related 
to the spin polarization, which suppresses the spin fluc¬ 
tuations. In the field along the [110] direction, neutron 
scattering has revealed a field-induced transition at 0.5 
T42 With field lower than 0.5 T, there is significantly dif¬ 
fuse scattering intensity originated from the spin fluctua¬ 
tions. With increasing field, the diffuse scattering inten¬ 
sity is weakened and well-defined spin-wave excitations 
appear at 0.5 T, indicating a long-range ordered state 
(spin polarization) i 43 Since the spin waves are gapless at 
the critical field of 0.5 T, they can be well populated 
by the thermally excitations even at very low tempera¬ 
tures and strongly scatter phonons. This can explain the 
sharp 0.5T-“dip” of n(H) curves. Although there is no 
experimental investigations on this field-induced transi¬ 
tion with field along other directions, the 0.5T-“dip” of 
k(H) curves with H || [111] may have the same origin. 


H. High-field k(T) of Yb 2 Ti 2 0 7 

In most cases, the magnetic scattering can be sup¬ 
pressed by strong magnetic field and the k would be 
independent of field in high fields and display a T 3 de¬ 
pendence at low temperatures. However, Yb2Ti 2 0 7 was 
found to be exceptional. Figure 11 shows the k(T) with 
8 T along three different directions. It can be seen that 
these data do not exhibit the expected T 3 behavior even 
at T < 100 mK; instead, the data follows a T 2 5 de¬ 
pendence at temperatures lower than 200-300 mK. The 
phonon mean free paths are calculated and are found to 
be smaller than the sample sizes, as shown in Fig. 11(d), 



FIG. 11: (Color online) (a-c) Temperature dependencies of 
the thermal conductivity of Yb 2 Ti 2 0 7 crystals in zero field 
and in 8 T along different directions. The dashed lines in¬ 
dicate that the low-temperature data in 8 T follow the T 2 '° 
dependence, (d) Temperature dependencies of the ratio of 
the phonon mean free path l to the averaged sample width 
W, calculated from the 8 T data. 

which means that there are still some kind of microscopic 
scattering on phonons and the boundary scattering limit 
is not achieved. This temperature dependence of n in 
8 T field is a peculiar phenomenon that is not easy to 
understand. Since high field can effectively suppress spin 
fluctuations, there should be some other factor damping 
the phonon transport. 

IV. DISCUSSION AND SUMMARY 

Firstly, there is no strong evidence for the large 
heat transport of monopoles in the spin-ice compound 
Dy 2 Ti 2 0 7 . Actually, most of the experimental data, in¬ 
cluding the zero-field k(T), the k(H) with different field 
directions, the high-field behaviors, and even the relax¬ 
ation phenomenon, are more or less incompatible with 
the supposition of large monopole transport in zero field. 
The quantitative analysis using some phenomenological 
theory gives a qualitative description of the low-field de¬ 
crease of k for H || Jh , but it fails to explain the nearly 
isotropic behavior with H || [111] for k along the field 
and perpendicular to it 42 Furthermore, the calculations 
show that very high speeds (> 10 4 m/s) of monopole 
propagating are required to describe the experimental 
data. This is unlikely since the monopole excitations are 
dispersionless and in principle they can only transport 
diffusively with a small averaged velocity. In the case of 
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Yb 2 Ti 2 C> 7 , the kink-like feature of k(T) at 200 mK is a 
clear signature of considerable monopole transport at the 
quantum spin-ice state (T > 200 mK). Since the quan¬ 
tum monopoles are dispersive, it is possible they have a 
highly mobile characteristic. 

Secondly, judging from the k(T) data and the esti¬ 
mation of phonon mean free path, Dy 2 Ti 207 does not 
exhibit strong phonon scattering in zero field, while 
Yb 2 Ti 2 07 has rather strong scattering effect by the spin 
fluctuations. The high-field behaviors also demonstrate 
this difference. The main reason is that the monopole 
excitations are difficult to be thermally excited at very 
low temperatures due to the sizeable energy barrier (~ 
4.35 K). 

Thirdly, the field dependencies of n are much more 
complicated in Dy 2 Ti 207 . One phenomenon is that 
the low-field transitions have strongest effect and cause 
step-like decreases of k. Since both the qualitative and 
quantitative analysis based on the assumption of large 
monopole heat transport cannot describe the experimen¬ 
tal results, some unknown phonon scattering processes 
are relevant. In Yb 2 Ti 2 07 , the field dependencies are 
mainly determined by the suppression of spin fluctua¬ 
tions when the spins are polarized, and the field-induced 
transitions can cause a sharp dip-like anomaly in k(H). 

Fourthly, for Dy 2 Ti 2 C> 7 , the k in low fields exhibit a 
strong relaxation effect and irreversibility. These phe¬ 
nomenon are not observed in the k measurements of 
Yb 2 Ti 2 C> 7 . This is mainly related to the extremely slow 
spin dynamics in the low-field states of Dy 2 Ti 2 C> 7 . 

Fifthly, there are some peculiarity in the high-field heat 
transport of both materials. Dy 2 Ti 2 C >7 displays strong 
field dependencies of k in high field up to 14 T, which 
is quite anomalous considering the spins are polarized 
above ~ 2 TjiS Yb 2 Ti 2 C >7 shows field independent k in 
the polarized state. However, the T 25 behavior of high- 
field k(T) is not expectable if the magnetic scattering 
of phonons are completely suppressed. These indicate 
that besides the impacts associated with the low-field 
magnetic transitions and spin fluctuations, some other 
factors are also important. 

One factor can affect the phonon transport is the mag¬ 
netostriction and magnetoelastic coupling in these mate¬ 
rials. For example, at the transition from the kagome-ice 
phase to the saturated state with H || [111],— the ul¬ 
trasound measurement revealed sharp anomalies of the 
sound velocity and sound attenuation— Furthermore, 
there is a low-field phenomenon of Dy 2 Ti 207 that is 
not mentioned in the above discussions. For H || [100] 
or [110] and at very low temperatures, the weaker and 
broader step-like decreases of k at 12 T could not be 
related to any magnetic transitions and should have a 
different origin. In fact, all the earlier experimental re¬ 
sults of magnetization and neutron scattering have not 
found any magnetic transition or anomaly at 1-2 T for 


H || [100] or [110]. Instead, this phenomenon may be 
related to the magnetostriction, which exhibited a small 
anomaly at 2 T along the [100] direction— In this re¬ 
gard, the lower-field (< 1 T) transitions can affect the k 
in the same way since they show stronger effect on crys¬ 
tal latticed The high-field behavior of Yb 2 Ti 2 C >7 may 
be related to magnetoelastic coupling as well. In the 
spin-liquid Tb 2 Ti 207 , a magnetoelastic mode is formed 
by the hybridization of the first excited crystal-field level 
and the transverse acoustic phonons— As a result, the 
phonon transport of Tb 2 Ti 2 C >7 is so strongly damped 
that it behaves like a glassy stated In Yb 2 Ti 2 C> 7 , the 
coupling of the sound wave to quantum fluctuations 
has been investigated by the sound-velocity and sound- 
attenuation measurement. 83 Some anomalies in temper¬ 
ature and field-dependent sound velocity and attenua¬ 
tion were observed and were attributed to the first-order 
phase transition. Moreover, these measurements (down 
to 20 mK) have found that both the sound velocity and 
sound attenuation change continuously with field up to 
5 T, without showing any signature of saturation— This 
indicated that the magnetoelastic coupling is significant 
even in high fields. The structural distortions caused by 
the magnetoelastic coupling would prevent the phonons 
from transporting ballistically. 

In summary, the low-T thermal conductivity of spin-ice 
Dy 2 Ti 2 C >7 and quantum spin-ice Yb 2 Ti 2 C >7 are studied. 
The detailed temperature and field dependencies of k in¬ 
dicate that the mechanisms of heat transport in these ma¬ 
terials are rather complicated. The phonons are the main 
heat carriers and several factors including the magnetic- 
monopole excitations, the field-induced magnetic tran¬ 
sitions, the phonon scattering by spin fluctuations, and 
probably the structural distortions caused by the mag¬ 
netoelastic coupling are involved. One conclusion is that 
the monopoles in Dy 2 Ti 207 play a minor role, while the 
quantum monopoles in Yb 2 Ti 207 may make a sizable 
contribution in carrying heat. The comprehensive theo¬ 
ries for quantitatively describing the experimental results 
are called for and would be very useful for understanding 
the exotic magnetisms of these materials. Further exper¬ 
imental investigations on the magnetoelastic coupling in 
high fields are also indispensable. 
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